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Abstract 

We study a variant of the ferromagnetic Potts model, recently introduced by Tamura, 
Tanaka and Kawashima, consisting of a ferromagnetic interaction among q "visible" colours 
along with the presence of r non- interacting "invisible" colours. We introduce a random-cluster 
representation for the model, for which we prove the existence of a first-order transition for any 
g > 0, as long as r is large enough. When q > 1, the low-temperature regime displays a g-fold 
symmetry breaking. The proof involves a Pirogov-Sinai analysis applied to this random-cluster 
representation of the model. 
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1 Introduction 



Recently, in a series of papers [20l [2TJ [22] , Tamura, Tanaka and Kawashima introduced a variant of 
the ferromagnetic Potts model to study the relation between symmetry breaking and the order of 
the phase transition. The model consists of a ferromagnetic Potts interaction taking place between 
q "visible" colours along with the presence of r "invisible" colours without any interaction. They 
observed, through numerical simulations, that in two dimensions with g = 2,3,4 and r large, the 
model undergoes a first-order phase transition with g-fold symmetry breaking. This is in contrast 
with the ordinary two-dimensional g-colour Potts model with q = 2,3,4, in which the transition 
accompanied by the g-fold symmetry breaking is known to be of second order [I]. This provides 
a simple example that a q-fold symmetry breaking in two dimensions does not universally identify 
the order of the transition. 

The transition in this model (as well as in the standard g-colour Potts model) occurs from an 
ordered state (which has a favoured direction among q possibilities) to a disordered state (which 
has no favoured direction) when the temperature is increased. In the standard g-colour Potts model 
with q small, this transition is of second order (no latent heat at the transition point) whereas in 
the Tamura- Tanaka-Kawashima version of the Potts model, for the same values of q but r chosen 
sufficiently large, the transition is of first order (the system absorbs heat during the transition, 
without changing its temperature). 

For the standard g-colour Potts model, when q is large enough, there is a variety of difi^erent 
rigorous proofs that the transition is of first order [Ml |3l [171 [16] (see also [HI [10]). As announced 
in an earlier communication [5], in the present paper, we prove, by minor adaptations of the 
proofs in [El [15] , that when q + r is large enough, the Potts model with q visible colours and r 
invisible colours undergoes a first-order phase transition. The proof is based on an application of 
the Pirogov-Sinai method to a random-cluster representation of the model. 

The phase transition could be better understood if one thinks of a state of the system as a 
possible resolution of the confiict between order and disorder. The conflict should be resolved 
locally and in every region. In an ordered region, the neighbouring sites tend to take the same 
colour so as to minimize the energy, while in a disordered region, the neighbouring sites take their 
colours independently to maximize the entropy. To establish the resolution of the order-disorder 
confiict, one needs to take into account the disturbance present at the interface between ordered 
and disordered regions (the contours). 

For the standard q-colour Potts model, the order-disorder confiict is niftily depicted in the 
Fortuin-Kasteleyn (a.k.a. random-cluster) representation of the model [5] [TUl [HI [12] • In this repre- 
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sentation, order is associated with the presence of bonds between neighbouring sites and disorder 
with the absence of bonds. In the same spirit, we introduce a variant of the Fortuin-Kasteleyn 
representation for the Potts model with invisible colours. The advantage of this new formulation 
is that it admits a neat definition of the interface between ordered and disordered regions. Now, 
having two reference configurations describing complete order and complete disorder — the one 
with every bond present, and the one with every bond absent — as in [16) . we can apply the 
Pirogov-Sinai method [l5l fT9l [23 l [2] . 

In Section [31 we describe the model and recall the formulation of first-order phase transition 
in the Gibbsian setup. Section [3] is dedicated to the introduction of a variant of the random- 
cluster model and its connection with the Potts model with invisible colours. In Section [4. 11 formal 
definitions for contours are provided, and it is shown how to rewrite the partition functions of the 
model in terms of contours. These contour representations are then reduced, in Section [4. 2[ to two 
abstract contour models, on which standard techniques can be applied. In Section [5l starting from 
the two contour models, we obtain two approximations for the free energy of the Potts model with 
invisible colours. If q+r is large, each of these approximations turns out to be accurate in an interval 
of temperatures, one whenever order prevails and the other when disorder is dominant. The two 
intervals exhaust all the temperatures and have a unique common point, which is the transition 
point of the system. Finally, the above two approximations are used in Section [5] to prove a 
first-order transition at the transition point. The occurrence of the symmetry breaking at the same 
transition point then follows, using standard properties of the random-cluster representation, which 
are reviewed in Appendix IA.3I 

Acknowledgments. G.I. and S.T. thank NWO for support. We thank Roberto Fernandez for 
helpful conversations. 



Let L denote the two-dimensional square lattice, which we think of as a graph (§,B), where S 
denotes the set of sites (identified by Z^) and B the set of nearest neighbour bonds. In the {q, r)- 
Potts model, each site i g § is in one of {q + r) colours 1,2, ... ,q,q + 1, ... ,q + r. Therefore, a 
configuration of the model is an assignment of values from the set {1, 2, . . . ,q,q + 1, . . . ,q + r} to 
the sites in S. The {q, r)-Potts model [20l [21] [22] is expressed by the formal Hamiltonian 



where S{ai = aj < g) is 1 if Ui = crj < q and otherwise. Each pair of neighbouring sites that have 
the same colour a < q contributes with energy —1, while sites with colours a > q do not contribute 
to the energy. The first q colours are hence called the visible colours, and the rest the invisible 
colours. If there are no invisible colours (i.e, if r = 0), the model reduces to the ordinary Potts 
model with q colours. As in the ordinary g-colour Potts model, the {q, r)-Potts model has precisely 
q periodic ground-state configurations, in which every site has the same visible colour. 

Following the usual approach, we describe the system in thermal equilibrium via probability 
distributions on the space of all possible configurations of the model. The Boltzmann distribution 



2 Potts Model with Invisible Colours 



2.1 The Model 




(1) 
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on a finite volume A C L with boundary condition lo at inverse temperature (3 is defined by 

where H\{aA'-^\D) consists of a finite number of terms in the formal Hamiltonian ^ corresponding 
to the energy of a a and its interaction with the boundary condition w. Namely, 



which can be decomposed as a sum 



{ij}6B 
i e S(A) or j e S(A) 



Ha {aALo^c ) - i/^ (^A ) + ^A™ (^ac^ac ) , (4) 

where H'j^*^{(ta) involves the interaction terms within A, and i?^°"'^'*(CTAaj^c ) represents the terms 
corresponding to the interaction of A with its boundary. The factor Z'^ (A) is a normalizing constant 
— the partition function — making /it^ ^ a probability distribution. More specifically, the partition 
function of volume A with boundary condition uj is given by 

Z^{A) = ^e-'3^A"("^)-^^A™"'("^"Ac) . (5) 
If we ignore the boundary term, we obtain the free-boundary partition function of volume A: 

0"A 

This is the normalizing factor for the free-boundary Boltzmann distribution on A. A Gibbs measure 
on the space of all configurations of the infinite lattice system, at inverse temperature /?, is a proba- 
bility measure /i whose conditional probabilities for every finite volume A, given the configuration uj 
outside A, are given by the Boltzmann distribution ^. More specifically, 

fi{A and B)= [ ^i'^^ ^^(A) fiiduj) , (7) 

JB 

for every event A not depending on the colours of the sites outside A and every event B not 
depending on the colours of the sites in A (see e.g. [7])- It follows from a compactness argument 
that such measures exist at every temperature. However, when the temperature is sufficiently low, 
it is possible to have several distinct Gibbs measures. The multiplicity of Gibbs measures is then 
interpreted as the possibility of co-existence of distinguishable phases of the physical system (in 
this case, the possibility of spontaneous magnetization in q different directions). We refer to [7] for 
details. 

One way to obtain Gibbs measures consists of taking the thermodynamic limit of the Boltzmann 
distribution with or without a fixed boundary condition. For a visible colour k, let w'^ denote the 
configuration of the lattice in which every site has colour k. Let /i^ denote a Gibbs measure obtained 
by taking a weak limit of finite-volume Boltzmann distributions with boundary condition w*^ at 
inverse temperature /3, when the finite-volume grows to the whole lattice. Similarly, we obtain a 
Gibbs measure by taking a weak limit of free-boundary Boltzmann distributions. 



4 



For every n > 0, let A„ denote the {2n + 1) x (2n + 1) central square in the lattice, which we see 
as the subgraph of the lattice induced by the sites in [— n, n]'^. The pressure of the model is defined 

by 

= lim ^-^logZ^(A„) , (8) 

in which S{An) denotes the set of sites in A„. The function —jf{P) is the free energy per site. 
The limit exists and is independent of the boundary condition oj (see e.g. [131 [HI)- We would also 
get the same limit as in ([8]) if we used the free boundary partition function. The particular choice 
of volumes used above is not crucial, and can be replaced by any sequence satisfying the van Hove 
property (see US]). 



2.2 First-Order Phase Transition 

A first-order phase transition in temperature is characterized by the presence of latent heat at the 
transition point 4 . This means that at the transition point, the system absorbs or gives out heat 
without a change in temperature. The presence of latent heat, therefore, corresponds to a jump in 
the internal energy. 

In the Gibbsian setup, the state of a system in thermal equilibrium is represented by a Gibbs 
measure (see [7[). If the Gibbs measure is translation-invariant, the internal energy density of 
the system is described by the expected value of energy per site. The presence of latent heat at 
a temperature means that the limits of the internal energy from above and below the transition 
temperature are different. This implies, by continuity, the existence of two translation-invariant 
Gibbs measures at that temperature having different internal energy. 

If the pressure function /(/?) is differentiable at a point /3, its derivative at (3 coincides with 
the internal energy with respect to every translation- invariant Gibbs measure at /3 (see [13|, I18j). 
(This, however, does not rule out the possibility of the existence of several translation-invariant 
Gibbs measures at /3.) If, on the other hand, the pressure function /(/3) is non-differentiable at a 
point /3, its left and right derivatives at /3 (which exist due to convexity) are different and coincide 
with the internal energy with respect to two different translation- invariant Gibbs measures at /3. 
The difference between these two derivatives corresponds to a latent heat at /3, implying that the 
system undergoes a first-order phase transition at /3. 

In this paper, we show the existence of a first-order phase transition in the {q, r)-Potts model 
by proving that the pressure function /(/3) has a unique non-differentiable point (3^- Above (3^, 
the system admits q "ordered" translation-invariant Gibbs measures. Each ordered measure can 
be thought of as a perturbation of one of the q ground-state configurations, in the sense that with 
probability 1, the configuration of the model consists of a unique infinite sea of one of the visible 
colours with finite islands of disturbance. Below /3c, the system has a "disordered" translation- 
invariant Gibbs measure, which can be seen as a perturbation of the uniform Bernoulli measure: 
there is a unique infinite sea of independent colours with finite islands of disturbance. At /3c, the 
q ordered measures co-exist with the disordered one. 
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3 Biased Random- Cluster Representation 



In analogy with the standard Potts model [6l[T0l[8], it is possible to rewrite the partition function for 
the (g, r)-Potts model in terms of the partition function for a variant of the random-cluster model. 
While the former is a model defined on sites, the latter will be a model defined on bonds. The 
random-cluster representation of the Potts model allows for an elegant formulation of the intuitive 
concepts of "order" and "disorder" : the presence of a bond in the random-cluster representation is 
interpreted as "order" , while the absence of a bond as "disorder" [H] . 

Although for the purpose of our problem, it suffices to present the connection for squares A„ in 
the lattice, we elucidate the connection for an arbitrary finite graph, where there is no boundary 
condition. Later, we explain how the boundary conditions affect this connection. 

Let G = {S,B) be a finite graph. The r-biased random- cluster model on G is given by a 
probability distribution on the sets X Q B. The distribution has three parameters 0<p< l,g>0 
and r > 0, and is defined by 



p.q.r \ 



heB 



in which kq (S, X) denotes the number of isolated sites of the graph {S, X) and ki {S, X) the number 
of non-singleton connected components of {S, X) and Z^^ ri'^) the partition function. Notice that 
for r = 0, the model reduces to the standard random-cluster model, in which both singleton and non- 
singleton connected components have weight q. For r > 0, the above model induces a bias towards 
singleton connected components. Namely, the singleton connected components have weight {q + r) 
whereas the non-singleton connected components have weight q. 

Let us now see how the (g, r)-Potts model is related to the r-biased random-cluster model. This is 
a mere generalization of the standard relation between the Potts and random-cluster models (TUl [5] . 
Let O be the set of {q, r)-Potts configurations on G. The partition function of this model can be 
rewritten as 

Zp{G) = e''^{'-}E«^("'="^-«) (10) 

= E n e^'^^'^-'-"^ (11) 
cren {i,j}eB 

= E n [^ + Si'^^ = '^J<1)(^^-^)] (12) 

o-eo {i,j}eB 

= EE(^'-1)"'' n Sia,=a,<q) (13) 
= E E '^(^'^) ' (14) 



o-eo xcB 



where 



7:{a,X)=S^\ n [S{{^,J}eX)S{a.^a,<q){l-e-^)+Si{^,J}(^X)e-^] . (15) 

{i,i}6S 
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The latter expression can be seen as a coupling of the (g, r)-Potts distribution on 51 and a probability 
distribution on the space {0, 1}-^. The marginal of this coupling on {0, 1}^ is simply the r-biased 
random-cluster distribution 0p^,g,r with pp — 1 — e~^. In particular, the weight Tr{a,X) can also 
be expressed as 

^(a,X)=e^l^l-lj.^(a,X). [] [pf, S{{i, j} e X) + {1 ~ p^) S{{t, j} X)] , (16) 

where 

Fr ^ {(a, X):a,^ a, < q for aU {i, j} e X} . (17) 

The effect of the bias in the r-biased random-cluster model reduces to an increase in the number 
of compatible configurations with a given X, which is driven by a larger number of choices for the 
colour of those sites constituting the singleton connected components. In short, for each X C B, 
we have 

J2 IfMX) ^q^"''^'^^\q + r)^°^^'^^ . (18) 

The above coupling could be interpreted in either of the following ways [lOl [8] : 

I. We first sample a according to the {q, r)-Potts distribution. Then, we choose the elements of 
X from B, randomly and independently, as follows: for each bond £ B with = aj, 

we put {i, j} in X with probability pp; for each bond € B with 7^ aj, we do not put 

{i,j} in X. 

II. We first sample X according to the r-biased random-cluster distribution (j)pf,^q^r- Then, for 
each non-singleton connected component of (5, X), we pick a random colour uniformly among 
the visible colours, and colour every site in the component with that colour. Last, for every 
isolated site in (5, X), we choose a random colour uniformly among all the possible colours. 
(The choices of colours ought to be independent of each other.) 

We can now use to obtain 

Zf,{G)^eP\^\Z^l^^,{<G) (19) 

with pi3 = 1 ~ . 

For finite subgraphs of the infinite lattice, we will be using only two types of partition functions 
for the (q, r)-Potts model, namely the one with free boundary and the ones with homogenous 
boundary conditions. In the following, we see how the above two types of boundary conditions 
translate into the so-called disordered and ordered boundary conditions for the r-biased random- 
cluster model. Although, setting G = A„, equation ((T^ already provides a relation between the 
free-boundary partition functions of the two models, we will work with a slightly different relation, 
connecting the free-boundary partition function of the (q, r)-Potts model to a partition function for 
the r-biased random-cluster model that involves a boundary condition. This new relation will turn 
out to be more convenient in the sequel. 

The free-boundary partition function for the (g, r)-Potts model can be written as 
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A„+i 



(a) 



(b) 



Figure 1: (a) A configuration in A'^^^™'^. (b) A configuration in X™^^_^. 



wliere Zp^'^f°'^'^ {An+i) is tlie partition function witlr disordered boundary condition for tfie r-biased 
random-ciuster modei. Tlie latter is defined by 

^p"r''(A„+i) = J2 pfl(l-p^)l^(^"+^)\^l(g + r)''"(^(^"+^)'^)q''^(«(^"+^)^^) , (21) 

wliere 

-^A^"'' ^{XC B(A„+i) : X n (B(A„+i) \ B(A„)) - 0} ■ (22) 
Similarly, for the boundary condition uj'' we get 

zf{A„) = q-' ■ (e^ - i)-|iJ(A.+AA.)l . e/^|B(A„+.)l . ^RC.o.^(A„^^) , (23) 

wliere Zp^'°'l^{An+i) is the partition function with ordered boundary condition for the r-biased 
random-cluster model, which is defined by 

where 

-^aI+i = C B(A„+i) : X D B(A„+i \ A„)} . (25) 

By A„+i \ A„ we mean the graph obtained from A„_|_i by removing all the sites in A„ and the bonds 
attached to them. Let us remark that although mathematically X^^^™'^ is simply the collection 
of all subsets X C i?(A„), we wrote it as above to emphasize that the elements of A'^'*'™'* are 
configurations of i?(A„+i). See Figure [T]for typical examples of elements in ;fdisord <%a"^'^j- 

In the following section, we will extend the definition of Z^^-^f"'^'^ and Z^'^°^^ to arbitrary 
subgraphs of the lattice. 
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Using the above relationships, we obtain that the pressure of the (g,r)-Potts model can be 
written as 



We call the limit 



|B(A„)| / \ n^oo |B(A„)| 



/'^^(^) = hm TBTm ^ 1™ IP,: M (27) 



n^oo |S(A„)| n^oo |S(A„)| 

the pressure (per bond) of the r-biased random-cluster representation. Note that the singularities 
of the {q, r)-Potts pressure function can be detected by studying the pressure function f^'~{f3). 
One advantage of this random-cluster representation is that it has a more transparent expression 
in terms of "contours", which helps us study the function f^'~{f3)- 



4 Reduction to Contour Model 

4.1 Contour Representation 

Any configuration of the r-biased random-cluster model in a volume is a subset X of bonds in 
the volume. We interpret each bond in X as "ordered" and each bond outside X as "disordered" . 
Any configuration X can then be seen as clusters of ordered and disordered bonds. Whether an 
equilibrium state is ordered or disordered can be seen as the result of a competition between ordered 
and disordered regions. The selection criterion for this competition is "energy" . The term "energy" 
refers to an abstract notion of energy for the r-biased random-cluster model, which in analogy with 
the Boltzmann distribution, corresponds to minus logarithm of probability. 

Let us define the "energy" of an ordered bond as the "energy" per bond of the fully ordered 
configuration; that is, 

-loe \a(l - e-^'V^f^")!] 
e(B) ^ hm ^ = - log(l - e-'') . (28) 

n^oo \B[An)\ 

Similarly, define the "energy" of a disordered bond as the "energy" per bond of the fully disordered 
configuration: 

e(0) ^ hm = - b-'V^] • (29) 

n^oo |^(A„)| 

The "energy" of the ordered and the disordered regions can now be expressed as |i?°| • e(B) and 
\R'^\ -6(0), respectively, where |i?°| and \R'^\ denote the size of the ordered and disordered regions. 
The "energy" of X, in turn, can be written in terms of the "energy" of ordered and disordered 
regions plus a correction term due to the effect of the boundaries separating them. If the effect 
of these boundaries is negligible (which will turn out to be the case whenever g -I- r is large), the 
selection criterion for the competition between order and disorder boils down to determining which 
of e(B) and e(0) is minimal. This is the starting point of the Pirogov-Sinai approach to study 
phase transitions (see e.g. [H]). 
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The presence of the correction term at the boundaries can be explained as foUows: In the 
probabihty weight of a configuration X, each isolated site contributes with a factor {q + r). To 
express the "energy" of the disordered regions purely in terms of bonds, we evenly distribute the 
contribution of the isolated sites among the 4 incident bonds. Doing so, every disorder bond acquires 
zero, one or two "energy" -shares, depending on the number of isolated sites it is incident to. Since 
in the fully disordered configuration there is no ordered region, every bond is incident to precisely 
two isolated sites and receives two "energy" -shares, leading to the factor {q + r)^ in the expression 
of e(0). In an arbitrary configuration, however, the disordered bonds on the borderline between the 
ordered and disordered regions, receive one or no "energy" -share, hence the need for a correction 
term. 

It is possible to define a suitable notion of boundary between ordered and disordered regions, 
so that each configuration X is uniquely identified by its boundary (see below). We could then 
rewrite the partition functions as sums running over "admissible" boundaries, that is, those cor- 
responding to configurations of bonds. Each admissible boundary is split into "primary" objects 
termed contours whose "energy" add up to the corresponding boundary effect. 

In the following, we specify rigorously the above heuristic notions of "boundary" and "contours" . 
We define the boundary of a configuration X C B as the set 

dX = {{i, b)eSxM: i^b and ie S{X) and b i X} , 

where i b means site i and bond b are incident, and S{X) is the set of sites incident to bonds in X . 
The set dX uniquely determines X. We say that two bonds b and b' in the lattice are co-adjacent 
if they belong to the same unit square. More intuitively, co-adjacency is equivalent to adjacency 
in the dual lattice. A set of bonds X is co-connected if for every two bonds b,b' e X, there is a 
sequence 6 = 61, 62, = ^' of bonds in X such that bi and fe^+i are co-adjacent. A contour is 

a set 7 C S X B such that 

i) the set of bonds appearing in 7, denoted by ^(7), is co-connected, and 

ii) there exists a configuration X such that {S{X),X) is connected and 7 = dX . 

We shall denote by F the set of all finite contours in L. If 7 is a contour, then removing the bonds 
5(7) breaks the lattice L into connected subgraphs. If 7 is finite, the graph L\ B{'j) has a unique 
infinite connected component, which we call the exterior of 7 and denote by ext7. The subgraph 
L \ 5(7) \ ext7 (which could be empty or disconnected) is called the interior of 7 and is denoted 
by int7. By V{'j) we will mean the union of int7 and the subgraph induced by i?(7)0 Let 7 be 
a finite contour. The configuration X such that {S{X),X) is connected and 7 — dX (which exists 
by definition) is either finite or co-finite. If X is finite, we call 7 a disorder contour, and if X is 
co-finite, we call 7 an order contour. Note that, if 7 is a disorder contour, all the sites appearing 
in 7 are in the interior of 7, whereas if 7 is an order contour, all the sites appearing in 7 are in the 
exterior of 7. As a result, we can safely represent a finite contour 7 by the pair (5(7), cc) where x 
is a label specifying the type of the contour (disorder or order) . This also means that the set of all 
finite contours T can be partitioned into two subsets: the set of disorder contours, which we denote 
by F'^, and the set of order contours, which we denote by F°. 

Two contours are said to be mutually compatible if they are disjoint (as subsets of S x B). Let 
us emphasize that two mutually compatible contours are allowed to share either sites or bonds, but 
not pairs. 

'^By the subgraph induced by a set of bonds we mean the graph obtained by those bonds and their endpoints. 
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Figure 2: (a) A disorder contour and its defining configuration, (b) An order contour and its 
defining configuration, (a') and (b') Geometric illustrations of (a) and (b). 
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If X C B is an arbitrary configuration, there could be several ways to partition its boundary dX 
into mutually compatible contours. One way to construct such decomposition in an unambiguous 
way is as follows; first, we partition {S{X)^ X) into its maximal connected components {S{Ci), Ci). 
Then dC'i form a partitioning of dX . Now, the maximal co-connected components of every Ci are 
contours that we identify as the contours of X. 

The above decomposition allows us to think of dX as a family of mutually compatible contours, 
which we call the contour family of X . Let us recall that the contour family of a configuration 
X uniquely determines X . However, note that not every family of mutually compatible contours 
corresponds to a configuration. In particular, in a contour family of a configuration X, between 
every two nested finite contours of the same type, there necessarily lies a contour of the other type. 
This requirement induces a long-range constraint among contours, which raises some difficulties in 
dealing with the contours. We will see later how to get rid of such a constraint. Let us call a family 
d of contours admissible if it is the contour family of a configuration X C M. We shall denote by 
A the set of all admissible contour families. A contour 7 in a mutually compatible family d of 
contours is said to be external if it is not in the interior of any other contour in d. Note that if d is 
an admissible contour family with no infinite contours, all the external contours in d are necessarily 
of the same type. 

Having formalized the notions of boundary and contours, we can now express the weight of 
a configuration of the r-biased random-cluster model in terms of the "energy" of its ordered and 
disordered regions and the correction term due to the contours separating them. The one-to-one 
correspondence between the configurations and the admissible families of contours allows us to 
write the partition functions as a sum over contour families. The ordered/disordered boundary 
conditions on the configurations translate into the constraints for the corresponding contour family 
that the outermost contours in the volume be of the order/disorder type. 

Let A be a volume in the lattice, by which, from now on, we shall mean a finite subgraph of L 
without "holes" . More precisely, we assume that if we remove the subgraph A from L, the remaining 
subgraph is connected. Let us denote by AJ^"^"'''^ the set of all admissible contour families whose 
contours are in A (i.e., their bonds are chosen from the bonds of A) and whose external contours are 
all of the disorder type. Similarly, let A""^ denote the set of admissible contour families in A whose 
external contours are all of the order type. The partition function for the r-biased random-cluster 
model in a volume A with disordered (resp., ordered) boundary conditions can be defined as 

<,fr'(A) ^{q + Yl e-I^A(^)l-W-|^A(a)h(0) -Q ^(^) 



(31) 



where ^(7) is the weight of a contour 7 and is given by 

/ A A 1(9 + '")^^''''', if 7 order, 

\qyq + r) 4^1^ it 7 disorder, 

and R\{d) and R%{d) denote, respectively, the sets of ordered and disordered bonds in A of the 
configuration corresponding to d. 

The above definitions are consistent with the definitions given in (f2T|) and ([24|) when A — A„+i 
is a square. Namely, for A = A„+i, if X is the corresponding configuration of a family d G A^'^°'''^, 
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the restriction of X to i?(A„-|_i) is an element of X^^^™'^. Conversely, every element of 
has a unique infinite- volume extension whose corresponding contour family is in A^"'°'^'^. A similar 
correspondence holds between '^a^'^^^ and A^''^^^. For the proof of the equivalence of the two 
definitions see Appendix I A. 2 1 

|aB(A)| 

We emphasize that the factors (q+r) J and q in front of the partition functions ((3T|) and p2l) 
do not contribute to the pressure function f^^{f3). in the thermodynamic limit, they are swallowed 
by the size of the volume. Hence, to avoid heavy notation — with all due apologies to the reader — 
we re-define the partition functions of the r-biased random-cluster model with disordered/ordered 
boundary conditions as 

^RC.disord^-^^) ^ Y Q-\Rlid)\-e(«)-\Ri{9)\-e{0) Y^pi^^^ ^ ^3^^ 

From now on, every time we talk about the partition function of the r-biased random-cluster model, 
we will be referring to the latter definitions. 

As was mentioned in the introduction, we would like to express the two partition functions in 
terms of two (standard) contour models. The purpose of this is to make use of the machinery 
available for contour models; namely, a result providing an estimate on the convergence of the 
corresponding free energy functions (Proposition [T]), and the Peierls estimate for the probability 
of the appearance of a contour. The main features of the contour models that are required in the 
above tools are (see [TS]) 

i) independence, and 

ii) damping. 

Unfortunately, the contours of the contour representation of the r-biased random-cluster partition 
functions are not independent (due to the long-range constraint). In the following section, we 
will see how to achieve the independence among contours, by rewriting the partition functions in 
terms of abstract contour models. As in the standard random-cluster model (see [IS]), we need two 
different such contour models, one for each of the two boundary conditions. 



4.2 Contour Models 

In this section, we want to resolve the issue of long-range constraints between contours. Recall that 
the admissibility condition requires the contours of a family to be alternating between disorder and 
order contours, and this imposes a long-range constraint between contours. For example, two nested 
contours of the disorder type (no matter how far from each other) are "aware" of the presence of 
an order contour separating them. As a result, if we remove a contour from an admissible family, 
the admissibility could be lost. 

In order to get rid of this constraint, we use two abstract contour models in which the contours 
are all of the same type and the admissibility condition is replaced by mere mutual compatibility. 
The weights of the contours in each of the abstract models will be chosen in such a way to guarantee 
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that the ensuing partition functions are equal (up to a factor) to each of the partition functions for 
the r-biased random-cluster model. 

A contour model is specified by a function x : F — R, assigning a weight xil) to each contour 
7 e r. The configurations of the model are families of mutually compatible (i.e., disjoint) contours 
in L. Let us denote the set of all such families by Ai, and the set of all elements of A4 whose 
contours are in a volume A by Ma- The partition function of the model in A is given by 

^(A|x)= E n^(^)- (36) 

deMA 769 

In the following lemma, we will see how to represent the partition functions of the r-biased random- 
cluster model with disordered and ordered boundary conditions, each in terms of of the partition 
function of a contour model, with a particular choice of the weight function. In fact, the contour 
model associated to the disordered boundary condition will not involve order contours. This is 
reflected by the fact that in this model each order contour has weight zero. Similarly, the contour 
model for the ordered boundary condition involves only order contours. 

To set the stage for the lemma, we rewrite the partition functions ^RC.disord^^^ ^RCordj^^Y) 
in a form resembling more the contour model partition function ^(A \ x)- That is, 

^RC.disord(^) ^ g^|S(A)|.e(0) ^ Y[p{j), (37) 

^RC.ord(^) ^ g-|i3(A)|.e(B) ^ [] ^(7) , (38) 

gg^ord^gg 

where 

~( _ jp'^^) ■ e-l^('"*'^)l-(^W-'=(^» , if 7 is disorder, 

~ \ph) ■ e-l^(^(^))l-('=(^)-^(»)) , if 7 is order. ^^^^ 

Let us recall that the set A^"^°'''^ (resp., A^'*^) does not contain only families of disorder (resp., order) 
contours, but all families compatible with the disordered (resp., ordered) boundary condition. To 
make the proof more transparent, let us define 

ag^ord ^ga 

so that 

^■RC.disord^^) ^ ^-|B(A)|e(0) . yd(^) ^ ^^2) 
^RC.ord(^) ^ ^-\BiA)\em . yo^^) _ (43^ 

Notice that the above contour representation for the partition functions ^RC.disordj^^-j ^RC.oidj-^-j 
lacks the condition of independence between compatible contours. 
The following lemma is similar to Lemma 1 of }15j . 
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Lemma 1. The partition functions for the r-biased random- cluster model on volume A with the 
disordered and ordered boundary conditions can be written as 



^KL.disord 
^RC.ord 



-|S(A)|e(0) I ^d^ ^ 

)ir(A|r), 

where the weights and ^° are defined by 



r(7) 



Pil) 
0, 



^RC.disord(iQ|-^) ' 



if 7 disorder, 
otherwise. 



(44) 
(45) 



(46) 



r(7) 



p(7)e 
.0, 



|B(7)|e(B): 



^RC.dii 



^RC.ord(ij^t^) 



if 7 is order, 
otherwise. 



Proof. The key step to prove the lemma is to write a recursion for the above partition functions by 
factoring the contribution of the interior of each external contour. Let us denote by f^isord 
of mutually compatible families of disorder contours whose elements are all external. (We include 
the empty family in Note that the elements of f^isord ^-^^ admissible and in A^"'°'''^. 

Moreover, for each admissible family d € A^"'°'''^, the sub-family of d consisting of its external 
contours is in Similarly, we denote by £^'^ the set of mutually compatible families of order 

contours whose elements are all external. The partition functions and Y° satisfy the following 
recursions: 



r°(A) 



E 



Similar recursions hold for the contour model partition functions | ^'^) and | ^°): 



1 766 



Note that, since every order contour is weighted by we can ignore in | ^'^) the families 
containing order contours, and similarly the disorder contours can be ignored in | ^°). 

We use induction on volume A to prove that Y'^{A) = iF(A | ^'^). Suppose that for every sub- 
volume A' C A we have y'i(A') = 3f {A' | Let 6* e £f''°''^. We want to show that the terms 
corresponding to 9 in the recursion formulas (|48l) and (|50|) for Y'^{A) and 3f{A \ ^'^) are equal. If 
is empty, the equality is trivial (we consider the product over an empty set to be 1). Otherwise, 



(47) 



(48) 



(49) 



(50) 



(51) 
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for every j G 9, we have int 7 C A, which iraphes 3f{int 7 | ^'^) — y^(hit 7). Usmg the definitions of 
p and ^'^ we obtain that 

n P(7) • r°(int7) = H eh) ■ Y%ntj) . (52) 



Therefore, y^(A) — ^(A|^'^). The starting point of the induction is when the only element of 
£disord ^Yie empty family 

The argument for Y°{A) = 3f{A \ is similar. □ 



Note that there is no complete correspondence between the configurations of the r-biased 
random-cluster model and the contour families of the corresponding abstract contour model. Nev- 
ertheless, the probability of appearance of a contour as an external contour is the same in both 
models. Let 0™*^ denote the probability distribution associated to Z^'--°^''^{A). We consider (j)'^'^ 
as a measure on the infinite-volume bond configurations X CM, which is concentrated on the set 
{X : dX e A"'^}. Likewise, 4>f''°'"^ will denote the measure corresponding to Z''^ '*'''°'''^(A), which 
is concentrated on the set {X -.dX e A^'''°'''^}. 

Corollary 1.1. Let A be a finite volume and 9 G E]^'^ a family of external mutually compatible 
order contours. Then, 



n.f=«r(7)^(int7ir) 
: ^extX = 9}= " ' ^''^ 



where 9cxt^ is the family of external contours of X . A similar statement holds for the probability 
of families of external mutually compatible disorder contours under . 



Proof. 



{X : 9e.tX = 9}= ^^^"^ (54) 



Yl.ee U7) •e-l^(^(-))l-(^(--^(»)))^:M] y°(int7) 



r°(A) 

n7e^e(7)^(int7ir) 

ir(A|e) 



□ 



7erA 76rA 

7 order 



(55) 
(56) 



The next corollary provides an estimate for the probability that a finite region of the lattice is 
surrounded by a contour (see e.g. [IQ]). For a finite set of sites A in the lattice, let Ta denote the 
set of all finite contours that have A in their interiors. 

Corollary 1.2. For every finite volume A and every finite set A C S'(A) we have 



A similar bound holds in the disordered case. 
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Proof. Taking into account the ordered boundary condition, we have that if A is surrounded by a 
contour in A, it is also surrounded by an external order contour in A, that is, 



{X ■.dxnTA^0} = {x : d,^tX n ^ 0} . (ss) 

By the previous corollary, we can bound the probability of a contour 7 appearing as an external 
contour by 

r(7)^(int7ie)Eee£r T (7)ir(int 7 | e) 

_ 697 



^(Aie) 



We obtain that 



(60) 



r(7)^(int7in^(A\int7ir) 

ir(A|e) 

< r(7) ■ (62) 

The last step follows from the fact that all the terms in the partition function iF(A | are non- 
negative, hence 



^(A I > J^(int 7 I $°)ir (A \ int 7 I C) ■ (63) 



r^^ix ■.dxnrA^0} = ^a^Ix d.^tX n ^ 0} (64) 

= J2 ^T{X : d^^tX 3 7} (65) 
76rA 

< E ^°(^) ■ (66) 

□ 

A standard argument using the positive correlation property of (pj^'^ (resp., cf)']^^"'^'^) can be used to 
show that the thermodynamic limit of 0^'^'^ (resp., <^^"'°'''^) exists and is unique (see Appendix IA.3p . 
The limit measure (j>°'^'^ (resp., ^^isord^ satisfies the same bound as in the above corollary. If the 
weights ^° (resp., ^'^) decay sufficiently fast, the latter bound implies that under (jf"^'^ (resp., 0'^"'°'''^), 
the configuration of the model almost surely consists of a unique infinite sea of order (resp., disorder) 
with finite islands of disorder (resp., order). By a "sea" of order (resp., disorder) in a random-cluster 
configuration we mean a connected component of present (resp., absent) bonds. 

Corollary 1.3. For every finite set A QS we have 

0-d{x : M n ^ 0} < J2 = E ^°(^) ■ (67) 

7 order 
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Furthermore, if the sum X^^gr^ C°(7) converges, we have 

^ord f 3 unique infinite sea of order \ _ 

with finite islands of disorder J ' 

A similar statement holds in the disordered case. 

Proof. As before, let A„ denote the (2n + 1) x (2n + 1) central square in the lattice. For every n 
let us define r^^A„ as the set of all contours in A„ having A in their interiors. Prom the previous 
corollary, we know that for every m, > n, the following bound holds: 

<PTl{x : dx n r^,A„ ^ 0} < Td) • (69) 

Since the event {X : dX fi F^.a^ ^ 0} is local, we obtain 

r^{X : dX n Fa,a„ ^0}<Y. • (70) 

due to weak convergence of (fP^^ to (fP^'^. Letting n ^ oo the first claim follows. 

^^^■y^VA ^"("T) converges, using a Borel-Cantelli argument, with probability 1, no infinite cas- 
cade of contours appears on the lattice. In particular, if we define 

S° — \ X ^ ^ unique infinite sea of order 1 

' with finite islands of disorder J ' 

S'^ — i X • ^ unique infinite sea of disorder 1 

y ' with finite islands of order J ' 

the latter implies that <p°"^{S° U S'^) = 1. We show that, in fact, (j)°'"^{S'^) = 0. 

Let j4 be a finite sot of sites in the lattice. For every X S one can find a volume A containing 
A such that d{X C\ B(A}) E AJ^'^"'*^ (i.e., the restriction of X to A is compatible with the disordered 
boundary condition). In particular, if we define 

C^j^^ ^ {X : 3 a finite volume A C A„ with S'(A) 2A&ndd{Xr\ B{Aj) e A^'*^"'^} , (73) 

we have Ca,Ai C Ca,A2 C • • • and S'^ C C^,a„- 

If m, n are integers with m > n, every configuration X that is compatible with the ordered 
boundary condition on A^ (i.e., dX € Aa^) and is in Ca,a„ necessarily has an order contour 
surrounding A. Therefore, by the previous corollary, we have 

^aI (Ca,aJ < rA^ {X:dXnTA^0}< J2 ?°(^) • (74) 

jerA 

Since Ca,a„ is a local event, by weak convergence of ^a^ to (p"'''^ we have 

-yeFA 
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Letting n — >■ oo, we obtain 

r'^iS^) < lini r"'{CA,Aj < V r(7) • (76) 

n— >C30 ^ — ' 

The latter holds for every finite set A C §, which by the convergence of the series, implies that 

5 Damping of Contour Weights 

One advantage of working with contour models is that when the contour weights are sufficiently 
"damped" (i.e., decay exponentially in the length with a sufficiently fast rate) the free energy exists 
and is bounded, and moreover, the error in the finite- volume approximations of the free energy is 
of the order of the size of the boundary of the finite volume. This is the message of the following 
well-known proposition (see e.g. Section 2 of [23^, or Proposition 2.3 of p[9]). 

Proposition 1. Let t > be sufficiently large, and suppose that the weight function % : T — !■ R o/ 
a contour model satisfies < xil) ^ e"'^''''! for every contour 7. Then, the limit 

= 1™ ir/a m log ir(A„ I x) (77) 

n^oo |i!(A„j| 

exists and satisfies g{x) < E7-s(7)90 ^(7) — In particular g{x) — > as r tends to infinity. 

Furthermore, there is a constant C = C (r) , such that C ^ as t goes to infinity, and for each 
finite volume A € Z'^ 

e9(x)|B(A)hC(r)|9A| < | < ^g(x)\B{A)\+C{T)\dA\ ^ ^.j,g^ 

where dA denotes the boundary of the volume A and can be defined as the set of bonds that are not 
in B{A) but are incident to A. 

The main purpose of this section is to identify conditions on the parameters (q + r) and /? under 
which the weights ^'^ and ^° are damped (i.e., satisfy the condition of the above proposition). We 
will see that when {q + r) is large, for any value of /3 > at least one of and ^° is damped, 
and moreover there exists a unique /? at which both ^'^ and ^° are damped. Let us remark that 
for sufficiently damped weights, the sum appearing in CoroUarv 11.31 converges, implying that the 
corresponding phase is stable. 

For T > large enough, let us introduce the truncated (i.e., artificially damped) weights 

^-o(,) [^(7), if^°(7)<e-^H, ^^^^ 
I , otherwise, 

and similarly for ^^^(7) (see e.g. 115 ). The term truncated refers to the suppression of all contours 
whose weights are not damped. If we replace the original weight ^° by the artificially damped one 
we obtain the following truncated partition function, which can be thought of as an approximation 
of the partition function of the r-biased random cluster model with ordered boundary condition: 
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Similarly, replacing ^'^ by f*^ leads to the truncated partition function for the r-biased random 
cluster model with disordered boundary condition: 

^RC.disord(^) ^ g-|S(A)|.e(0)^(^ | ^d) ^g^^ 

The advantage of introducing these truncated partition functions is that we can apply Proposition[TJ 
Note that if the original weights are "damped" (that is, ^"(7) < e^"^''''! or ^'^(7) < e^'^l'^'l), the 
corresponding truncated partition functions coincide with the original ones. 

From Proposition [1] we have the following bounds for the truncated partition functions: 



5(9(?°)-e(B))|S(A)|-C(r)|aA| < ^RC.ordj-^^ < ^{g{e)-<B))\B{A)\+C{r)\dA\ ^ ^^^^ 
(g(5'')-e(0))|S(A)|-C(r)|aA| < ^RC.disord ^^^^ < g(g(«'')-e(0) ) | B(A) |+C(t) |aA| _ ^g^^ 



The pressure functions associated to the truncated partition functions are 

^°(^) = 1™ n?7rTi°g^^^"''''(^") = -<^)+9{n , (84) 

f{/3) = hm ^-^logZR^-d--<i(A„) = -e(0) + ) . (85) 

The functions and f^{f3) are lower approximations of the pressure f^'~{f3) of the r-biased 

random-cluster representation. The next lemma states that in fact when {q + r) is large enough, 
the maximum of f° and coincides with f^'~. As we will see in the next section, for {q + r) 
large enough, the functions g{£,°) and g{£,'^) and their /3-derivatives are small, and therefore, the 
dominant terms of f° and are — e(B) and — e(0). This means that f^'~ is approximated by the 
maximum of the curves — e(B) and — e(0), which intersect at a unique value /?, with significantly 
different slopes. 

By the diameter of a contour 7, denoted by diam7, we shall mean the maximum lattice distance 
between two bonds in -6(7). The next lemma is parallel to Lemma 2 or [15^ or Theorem 3.1 of 

Lemma 2. Let {q + r) be sufficiently large. If < f° , then 
i) for every disorder contour 7 with diam7 < ^^^a we have 

< e-^l^l , (86) 

ii) for every order contour j we have 

Cil) < e-^l-^l . (87) 

A similar statement holds if f° < f^. 

Proof. We prove the two claims simultaneously by induction on diam7. Let K > and suppose 
that the claims hold for all (disorder/order) contours with diameter less than K. 
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Let 7 be a disorder contour with diameter K that satisfies diam 7 < fa^fd ■ Then, 

.d. N . N ^'^'^•°'''^(int7) 



^-RC.disord^jj^t^) 
^RC.ord(ij^^ ) 

-^(7) ^RC.diso.d(i,t7) ^''^ 

g/°-|B(int7)|+C(r)|aint7| 

= p(7)e(^°"-''''^l^('"*'^^l+^'^(^)l^'"*^l , (91) 

where in the second equahty we have used the induction hypothesis. Namely, every contour in int 7 
has diameter less than K, allowing us to replace the original partition functions with the truncated 
ones. Notice that 

• Pil) = qiq + ryi\'^\, 

• |i3(int7)| < i I7I • diam7, 

• if° - •diam7 < 1, and 

• |9int7| < I7I. 
Hence, we obtain that 

^''{1) < qe-(i log(g+r)-l-2C(r)).|7| _ (92) 

For (q + r) large enough (uniformly in 7) the righthand side is bounded by e^'^''*'! , hence the claim. 
Next, suppose that 7 is an order contour with diameter K . We need to show 

= oMel^^-'^I^W ^RC.disord(y(^)) 

? W Pl7je ^Rc.ord(int7) " ' ^^^> 

By the induction hypothesis, the partition function Z'^*" °'''^(int 7) is equal to the corresponding 
truncated partition function, which can be bounded using Proposition m As for Z^'^-'^'''°"^{V{j)), 
if we suppress all the contours that are "big" , we can get a similar bound using the induction 
hypothesis. 

To render the argument more transparent, we work with the partition functions Y'^[A) and 
r'i(A) (see go]) and dH])) for which we have 

^RC.disord^^) ^ g-|S(A)|.e(0)yd(^) ^ ^g^^ 

ZRC--<1(A) =e-l^(^)l-(»)y°(A) . (95) 

Let us call a disorder contour small if its diameter is less than or equal to jarrp- Otherwise, we 
call the contour 6117. 

As before, let us denote by £'disord ^j- mutually compatible families of disorder contours 

in A whose elements are external. Factoring the contribution of the interior of big external contours 
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we have the recursiorH 



F^(A)= J2 i"smaii(A\int0) nP(7')>^°(int7') (96) 
e big 

where iiit6l = |J^,ggint7', and >;iau(A) - ^(^Ksmaii)' in which the weight function ^^^^^ji is 
obtained from by replacing the weights of all big contours with 0, that is, 

10 if 7 big. 

The expression for the weight ^°(7) reads then 

e big 

The induction hypothesis and Proposition [T] tell us: 

• ^slall(^(7) \ intfi*) < e9(it,^n)-\B{V{'i)\\ryt0)\+C{Ty\d{y{'y)\\Yyte)\^ 

• r°(inte') < e9(n-\B(inte)\+c(r)-\dinte\^ 

• r°(int7) > Q9ii°y\B{int-i)\-C{ry\dinti\ ^ 

• g{^°) < e~^/2 < 1. 

Moreover 

. \d (F(7) \ int 0)1 < \dVm + laint^l < 3 I7I + Eyee IVI, and 

• 1^(7)1 < l7l- 
Hence we have 

r(7) < P{l) ■ e(i+''^(^»-l''l E el^(^('^)\'"*'')l-['=W-'=('2')+9(«?maii)-9(r)] p(y) . e2C(r).|7'| ^gg^ 

e big 

= p(7) • e(i+4<^M)-l^l E el^(^(^)\'"*'')l-[(f(«'man)-9(«"))-(r-/")] Y[ p{i) ■ e2C^(^)-h'l . (lOO) 
gggdiBord ^/ge 
e big 

As we shall see shortly, if {q + r) is large enough, the sum appearing in the above expression can 
be bounded by e^'^^'^' l'*'! , so that 

r(7) < P(7) • e(i+4C^M)-l''le3C^(^>l^l (101) 

= e~(3 log(9+'-)-l-7C(T)).|7l (102) 



^Although A\int6 does not match our requirement for being a volume (i.e., not having holes), it does not cause 
any problem. In fact, since the contours in ^gf^jj^u are small, they cannot surround the holes in A \ int6, hence they 
do not distinguish the holes from the outside of A. 
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For large {q + r), the righthand side is bounded by e '^''''l, proving the claim. 
It remains to show that for (q + r) sufficiently large, 



Q\BiV{j)\inte)\-[{g{it^^^,}^g(i''))-{r-f'')] "Q ^(y) . g2C(r).|7'| < e3C(r).|7| _ (^33^ 

e big 

To show this, let us consider a contour model with weight function 

/^ fp(7') ■ e^'='(''' l'''l , if 7' big and disorder, 

P(7 ) = < '^^ ^ ' ^ (104) 

I , otherwise. 

Assuming that g{p) < f° - + .g(C'*) - 5(CLaii)' and (g + r) in such a way that ^(7') < e~'"l'^'l, 
we can use Proposition [1] to obtain 

Q\B{Vh)\int0)\-[{g(if^^„)-g{i''))-{r-f'')] "Q ^(-y^ . e2C(r).|7'| 

big 

< e-l^^^f'^)^"^*^)!-^^?) n ^(^')-«'^'"^''^'' 
e big 

= e-l^(^('^»l-9(?) ^ J]^ p(y) .e2C(^)-|7'| .e|sM'r')|-9(p) (loe) 
e big 

< ^-mvmigip) Y n ^(V) • e'^("HV| . ^(inty I p) . eC^M-b'l (107) 

e big 

= c-l^(^(^»l-9(^) • 3f{Vi-f) I p) (108) 

< e^^f^)-!'^! . (109) 

Finally, to see that_g(p) < /° - /'i + g(e'i) - 5(4maii)> note that ^(^"d) - 3(^1^11) > (due to 
the fact that ^^maii — ^^d by Proposition [1] 

gip)< E P(7') < e-^/(^°-^') , (110) 

7':S(7')30 
7' big 

using the fact that 7' is big only if |7'| > y^zjd- For r not too small we have < 

From the above lemma, we know that the pressure f^'~{f3) of the r-biased random-cluster rep- 
resentation coincides with maximum of the functions and /'^(/3), that is, 

/R^(/3)=max{r(/3),/'i(/3)} . (ill) 
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(a) 



(b) 



Figure 3: (a) The curves of f°{fi) and f'^{fi) are within a narrow margin above — e(B) and — e(0). 
(b) The slopes of and f^iP) are close to the slopes of — e(B) and — e(0). 



Recall that 



-e(0)+5(^"'i) 



(112) 
(113) 



If T is large, Proposition[l]says that g{S,°) and g{S,'^) are small, so that f^'~{P) can be nearly expressed 
in terms of the "energy" per bond of the fully ordered and fully disordered configurations. More 
precisely, if we define 



F(/3) = max{- 



-e(0)} 



we have 



0</R^(/3)-F(/3)<e-^/^ 



The two curves — e(B) and — e(0) (as functions of /3) intersect at a single point 

= log (1 + x/gT7) , 

above which F{(3) — — e(B) and below which F{13) = — e(0). Furthermore, these two curves have 
significantly different slopes, implying that F{f3) is not differentiable at /3c. What we are after is to 
infer that f^'~{[i) has a similar behavior. In other words, we would like to show that there exists 
a unique solution /3c for the equation f°{/3) = /'^(/3), at which f^'~{/3) is not differentiable, above 
which = f°{l3) and below which /'^^(/3) = f^iP). Note that condition (fUS]) guarantees that 

/RC(^) = jo^^) > jd(^) foj. /3 > ^c and f^'^{/3) = > f°{l3) for l3 < ^c- In fact, it states 

that f^'^iP) lives in a margin of width e""^/^ above i^(/3) (see Figure I^Ka)). To infer such a sharp 
transition we further need to give bounds for the derivatives of g{^°) and g{^'^) (see Figure [Sfb)). 
This is addressed in the following lemma, which is analogous to Theorem 3.3 of [2]. 



(114) 



(115) 



(116) 
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Lemma 3. We have 



7:S(V(7))30 



~r/2 



,-r/2 



7:5(int7)90 

Proof. For a finite volume A, if we denote by r(A) the set of all contours in A, we have 

^iog^(Ain = ,;^^ E ^n?"°w 



dp 



^(^l^°)ok7ea 7.a 

7#7 



^ E 

7er(A) 

: E 

7er(A) 



d^il) ir(int7|C°) • Jr(Anext7|e) 



d/3 
dp 



^(A|e) 



Recall that ^"(7) is equal to either or C°(7), hence the derivative of ^"(7) is bounded by the 
derivative of C°(7)- Using the definition of ^°(7) (Eq. (H7|) ') we have 



9/3 



a/3 



dp ZRC-ord(int-y) 



The derivative of the partition functions appearing on the righthand side can be bounded directly 
from the definitions (Eq. ^ and ^) by 



d 



7-RC.disord 



(^(7)) 



< 



zRC.ord 



leading to 



Therefore, 



dp 

dp ZRC-ord(in|-^) 

dp 



(int7) < - |B(int7)| 



de( 



7RC.ord 



dp 



(int 7) , 



< 



-15(^7))! 



de(M) Z' 



RC.disord 



(^(7)) 



dp ZRC.ord(in|;^) ■ 



<-\Bivm\-^-m 
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We can now write 
d 1 



a/3|5(A)| 



logJ^(A|n <- 



< 



< -- 



< -2 



9e(B) 


1 


a/3 


\S{A)\ 


9e(B) 


1 


dp 


\S{A)\ 




1 


dp 


l^(A)| 


de{M) 


1 


dp 


\S{A)\ 


de{M) 


\B{A)\ 


dp 


\S{A)\ 



E \B{vm-c{i) 
E E 

7er(A) 6es(y(7)) 



a/3 



E E ^"(7) 
6eB(A) 7er(A) 

B(y(7))96 

E E 

6eB(A) 7er 

B(y(7))96 

E 

76r 

S(l^(7))90 

E ^7). 



7er 

S(V(j))30 



Since r is large, the above sum converges and leads to a uniform bound (with respect to A) for 
^ ^ log 3f{A I ^°). Using the dominated convergence theorem the first claim follows. The proof 



of the other claim is similar. 



(128) 



(129) 



(130) 
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□ 



The Main Result 

Theorem. Let e > and {q + r) large enough. The two-dimensional {q,r) -Potts model undergoes 
a first-order transition in temperature with breaking of permutation symmetry: 

i) Above the transition temperature, the model has a unique Gibbs state /x^'''^'', which is "disor- 
dered". 

a) Below the transition tem,perature, there exist at least q different "ordered" Gibbs states /i^, /i^, . . . 

Hi) At the transition temperature, q "ordered" Gibbs states /x^, /U^, . . . , /i'^ coexist with a "disor- 
dered" Gibbs state iJ'^^'^. 



The "ordered" and "disordered" states can be distinguished by 



/({ct ■.ai^k})>l-e, 
li'^^'^ia : a, = k}) < e , 



for every visible k, 
for every k. 



for every site i in the lattice. 



Proof. In the previous sections we have introduced the main ingredients to prove the occurrence 
of a first-order transition. Below, we first put these ingredients together so as to obtain a recipe 



(134) 
(135) 
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for the proof. Afterwards, we shall see how these ingredients along with basic properties of the 
biased random-cluster model (see Appendix IA.3[) can be used to prove the symmetry breaking at 
the transition temperature. 

The first step was to reduce the partition function of the (g, r)-Potts model to the partition 
function of the r-biased random-cluster model. This was done for the free and homogeneous visible 
boundary conditions, which led to the disordered and ordered boundary conditions for the r-biased 
random-cluster model (see Eq. (PO]) and ([25]) ). By means of this we could rewrite the pressure /(/3) 
for the {q, r)-Potts model as 

/(/3) =2/3 + 2/RC(/3) , (136) 

where 

f {(3) = Irni 1 = Irni |^77-T| ■ (137) 

The second step consisted in re-expressing the partition functions Z^'~ °'^'^ and ^RC.disord -^^ terms 
of the partition functions of two abstract contour models \ and iF(- 1 ^'^) (Lemma[T]) so that 
we could write 

f^{(3) = -e(B) + g{C) = -e(0) + g(^'^) , (138) 

where g{^°) and g{^'^) are the pressure functions for the two contour models (Eq. ([77]) ). 

If the weight function % of a contour model is sufficiently "damped" (i.e., xil) ^ e""^'"*"' for r 
large enough), the corresponding pressure g{x) can be made arbitrarily small (Proposition [T]) . In 
order to exploit this result, in the third step, we truncated the weight functions ^° and so as to 
render them artificially damped (Eq. ([79l) \ We could then define two functions 

/°(/3) --e(B)+g(e°) , (139) 
/^(/3) =-e(0)-f5(e') , (140) 

which approximate f^'~{l3) from below, and which can be thought of (for sufficiently large r) as 
perturbations of the functions — e(B) and — e(0), respectively (see Figure [Sj a)). 

In the fourth step, we proved that, for {q + r) large (relative to r), the pressure f^'~'{f3) is the 
maximum of these two approximations. This was achieved by proving that whenever f° > f^, the 
weight function ^° is "naturally" damped (i.e., ^° = and vice versa (Lemma [2]). Therefore, 
f^'~{f3) can be closely approximated by the maximum between — e(B) and — e(0). Due to the 
continuity of the pressure functions, the latter implies that the curves /°(/3) and f^{f3) intersect. 

In the last step, we showed that for r sufficiently large, -§pf°{P) > '§pf'^{P) (Lemma|3]). There- 
fore, the two functions meet at a unique point /3c at which f^''{(3) is non-differentiable (Figure[21[b)). 
Hence, the {q, r)-Potts model undergoes a first-order phase transition at /3c. 

We now prove the breaking of permutation symmetry at the transition temperature. Note that 
if /3 > /3c, we have f°{(3) > /'^(/3) and therefore, in view of Lemma [U the weights ^° satisfy 
£ e^"^!^!. Since r was chosen large, for every finite set of sites A, the sum X]7erA '^°('^) 
converges. Hence, it follows from Corollarv 11.31 that 

ord f 3 unique infinite sea of order . ^ 

Pf> V with finite islands of disorder ' ' ^'^^^'^ 
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iisord 



As in the previous case, Corollarv 11.31 and Lemma [5] guarantee that for every e > 0, choosing q + r 
large enough, ^^°°({(t : tXi = k}) < e. 

It remains to show that for f5 < j3c, the (g, r)-Potts Gibbs measure is unique. 

As in the standard random-cluster model, there exists a critical value < Pc < 1 such that 

• for p < pc, almost surely with respect to 0™'^ and 0p"^°'''^, there is no infinite connected 
component of bonds (order does not "percolate"), whereas 

• for p > pc, the event that a given site is in an infinite connected component happens with 
positive probability under both 0™'^ and 4>p^^™'^. 

(See Appendix lA. 31 ) It follows that Pc = pp^- Namely, if pp < p^, then order does not percolate 
under ^p'^'^. Therefore, equation (|14ip does not hold, implying that (3 < /3c. Conversely, if > Pc, 



(142) 



(The subscript is added to emphasize the dependence on /3.) Similarly, if /3 < /3c, the sum 
J2-yerA ^'*('^) converges and thus 

/ disord f 3 uniquc infinite sea of disorder 
Pi^ y with finite islands of order 

For a visible colour fc, let //^ be, as introduced in Section[2l a weak limit of Boltzmann distribu- 
tions with homogeneous boundary conditions w'^. Likewise, let p,^^°° be a weak limit of free-boundary 
Boltzmann distributions. 

Similar to the finite-volume couplings, there exists a coupling of /i^ and (p°'^^ with the property 
that with probability 1, every site incident to an infinite connected component of present bonds is 
coloured with k (see Appendix lA. 31) . If /? > /3c, we know that almost surely the bond configuration 
consists of a unique sea of order with finite islands of disorder. In particular, the probability that 
a given site i takes a colour other than k is bounded by the probability that site i is surrounded by 
an order contour; that is, 

: a, ^ fc}) < :dXnr,^0}. (143) 

In this region of (3, Corollary 11.31 and Lemma [2] ensure that (l)°'^^{X : dX n ^ 0} can be 
made arbitrarily small by tuning t. Hence, for every e > 0, choosing q + r large enough, we have 
^^({cr : (7, = fc}) > 1 - e. 

The measures ^^^'"^ and (t>p^^™'^ can also be coupled, in such a way that, given a configuration 
of bonds, the colour of the isolated sites are chosen independently and uniformly among the q + r 
possibilities. Using this coupling, and conditioning on whether a given site i is isolated or not, we 
obtain 

Ai|™({cr : cr, = k}) < ^— + (/)^|'"°''^{A: : i not isolated in (§, X)} . (144) 

p q + r 

If /3 < /3c, the bond configuration almost surely consists of a unique sea of disorder with finite 
islands of order. Hence, the probability that site i is not isolated is bounded by the probability that 
site i is surrounded by a disorder contour; that is. 



{X : i not isolated in (§, X)} < (pf^°"^{X ■.^XnT^^ 0} . (145) 
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then order percolates with positive probabihty under (pp^^™^, refuting (|142p . Hence, we must have 
/3>/3c. 

On the other hand, for every /3 at which 

0°'^'^(3 an infinite connected component of bonds) = , (146) 

the measure is the only Gibbs measure for the (q, r)-Potts model (see Appendix lA.Sp . The 
latter condition is guaranteed whenever < Pc, which is equivalent to /3 < /3c- Thus the uniqueness 
of Gibbs measure for (3 < [3^ follows. □ 



7 Conclusion 

In this paper, we presented a proof that the two-dimensional Potts model with q visible colours 
and r invisible colours undergoes a first-order phase transition in temperature accompanied by a 
g-fold symmetry breaking, provided the number of invisible colours is large enough. On the other 
hand, for r = (no invisible colours), the model reduces to the standard g-colour Potts model, 
for which it is known that if q — 2, 3, 4, the transition in two dimensions is second-order. Tamura, 
Tanaka and Kawashima [2Ql [211 122] introduced the Potts model with r invisible colours as a simple 
two-dimensional example with short-range interactions in which, tuning the parameter r, the same 
symmetry breaking could accompany phase transitions of different orders. The impossibility to 
infer the order of the phase transition from the broken symmetry was already noticed in other 
examples, such as the two-dimensional 3-colour Kac- Potts model f9]. For this model, Gobron and 
Merola proved that a 3-fold symmetry breaking might be accompanied with either a first-order or 
a second-order phase transition, by changing the finite range of the interactions. 

The first-order phase transition in the {q, r)-Potts model occurs as long as g -I- r is large enough. 
In particular, even for small values of q (say, q — 1,2,3,4), the presence of many invisible colours 
assures a first-order transition. The argument is very similar to the one for the standard q-Potts 
model, in which q is required to be large |16[ 115) . The transition point is asymptotically (in 
q + r) given by /3c ~ ■^\og{q + r). For q + r large, the latent heat is approximately given by 
2 (^^M^ + Mp.'^ ^2+ ^y2=, which tends to 2 as q + r oo. 

The proof relies on a formulation of the Potts model with invisible colours in terms of a variant 
of the random-cluster model, which we named the biased random-cluster model. The difference 
between this new model and the original random-cluster model is that it weights singleton connected 
components differently from non-singleton connected components. Such a disparity allows one to 
increase the entropy by increasing the number of invisible colours, while keeping the number of 
ground states (i.e., the number of visible colours) unchanged. The random-cluster representation 
allows for a clear formulation of order and disorder: order is associated with the presence of bonds 
while disorder with the absence of bonds. This leads to a simple notion of contours describing the 
interface between order and disorder. Hence, the random-cluster representation lends itself to a 
Pirogov-Sinai analysis, which is used to prove the existence of a first-order phase transition. 

We remark that the above analysis extends to higher dimensions. 
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A Appendix 

A.l Derivation of the Biased Random-Cluster Representation 

To derive the relation (f20| . we start from (|19|) and write 

Zp{K)=eP\^(^--^\-Z^l^^,{K) (147) 
Js:cs(A„) 

^ (g + ^)-|S(A„+AA.)l ^ (e/3_l)l>-|(5 + r)«o(5(A„+i),V)^«i(5(A„+i),n (149) 

ycs(A„+i) 
yn(B(A,.+i)\B(A,.))=0 

= (g + r)-l^(^"+AA„)l . e/^|s(A„^0l . Z^^^^fr^CA^+i) . (150) 

To obtain the relation (j23p . we need to take the homogeneous boundary condition for the (q, r)- 
Potts model into account. Denoting the set of (q, r)-Potts configurations on A„ by i^A„, we start 
from the definition ([5]) and write 



/3 ^ 5{a, ^<7j<q) + /3 ^ 5((7, a;}- = fc) [ (151) 
{jj}e-B(A„) {ij}eB(A„+i) 



»eS(A„),i^S(A„) 



<^6f2A„ {i,i}GB(A„) {jj}eB(A„ + i) 

ieS'(A„), J^S(A„) 



Denoting 



5A„ ^ e B(A„+i) : z G 5(A„) and j ^ S{K)} , (153) 

we can expand the products to obtain 

zf{K)^ E E E (e'^- 1)1^^1+1^^! (154) 
<TeOA„ XiCB(A„) X2caA„ 

where 

= {(7 6 Q.A„ Ci = (Jj < q for all {i, j} G Xi and cTj = fc for all i G ^(An) n S(X2)} ■ (155) 

To impose the ordered boundary condition, we multiply and divide by (e^ — l)l^(^"+i\A„)| 
emulate the presence of the bonds in i?(A„+i \ A„). This gives 

(A„) = (e^-1)-I^(^"+AA„)I ^ ^ (e'^ - 1)1^1 e efc(X)) , (156) 

<tGOa„ XC_B(A„+i) 

XDB(A„+i\A„) 
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where a is the extension of ct to a configuration in riA„+i with di — k for i e ^(An+i \ A„) and 

= {(7 e f^A„+i ■ o-i ~ dj < q for all {«, j} G X and di—k for all i G 5'(A„+i \ A„)} . (157) 
Changing the order of the sums gives 

zf{K) = {eP (e^ - 1)1^1 le^WI . (158) 

XCB(A„ + i) 
XDB(A„+i\A„) 

Note that for X satisfying i?(A„+i \ A„) C X C B(A„), the size of Qk{X) is 

q-l . (5 + ^).o(S(A„+i),X)^Ki(S(A„+i),X) (^5g) 

we obtain 

(A„) =<7'l(e^-l)-|^(A" + AA.)l ^ .(^^^)«o(5(A„ + i),X)g«i(5(A„ + i),X) (^gp) 

XCS(A„ + i) 
XDS(A„ + i\A„) 

= • (e^ - i)-|s(A.+AA.)l . e0|s(A„+oi . ZpR^^;°;'i(A„+i) . (161) 

A. 2 Derivation of the Contour Representation 

We show that for A = A„+i, the definitions ([21]) and dHU (resp., ([Ml) and ([32])) agree. 
The weight of a configuration X C _B(A„+i) is 

pPf|(l_p^)|B(A„ + i)\X|(^^^).o(5(A„ + i),X)^.i(S(A„+i),X) _ (^^2^ 

For a configuration X in '-t'^^'^°J''^ or <-f^'^^ ^ , let 9x be the corresponding contour family in A^'^°'j^ 
or A^';^^^. More precisely, dx = dX if Xf^'^^'^ and dx^d[x\J B{Kn+if) if -^aI+i- 

Claim. For X e we have the relation 

2 |B(A„+i) \ X| = 4ko(5(A„+i), X)+Y, \i\ - |9B(A„+i)| , (163) 

and for X G '^a^'^^ we have 

2 |S(A„+i) \ X| = 4ko(^(A„+i), X) + ^ |7l ■ (164) 

Proof. We first decompose the set {(i, b) : b E B{An+i) \ X and z ^ 6} as 

{{i,b) -.be B{An+i) \ X and i ^ b and i G S'(X)} (165) 

U 

{{i, b) -.be B{An+i) \ X and i ~ & and i ^ S{X)} , (166) 
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and furthermore note that the latter set can be expressed as 

{{i, b):ie 5(A„+i) \ S{X) and i ~ b} (167) 

\ 

{(i, 6) : i e S'(A„+i) \ S{X) and i - b and b B{A„+i)} . (168) 

We have 

\{{i, b):be B(A„+i) \ X and z ^ &}| = 2 |B(A„+i) \ X\ , (169) 
|{(i, 6) : 6 e B{An+i) \ AT and i - & and i G S{X)}\ = ^1 1^' ' 

|{(i, b):ie S{An+i) \ S{X) and z ^ &}| = 4ko(5(A„+i), X) , (171) 
and the cardmahty of {{i, b) : i £ 5(A„+i) \ S{X) and i ^ b and b ^ i?(A„+i)} equals 



|aB(A„+i)| , if X e A-^-^d^ 
,0, ifXeA-^:;^^. 



(172) 



□ 



Using the relations (|163p and (I164p the weight of X takes the form 

(g + r)3|9S(A„+i)|e-e(B)-|X|g-e(0).|B(A„+i)\X|(-^_^^-)-iE^6O^l7l^Ki(S(A„+i),X) ^^^3^ 

if X e A'f and 

g-e(B).|X|g-e(0).|S(A„ + i)\X|^^_^^^-iE-,ea^l7l^«i(S(A„ + i),X) ^^^^^ 

if X e '^a'"^j- If A^ G every non-singleton connected component in (5(A„+i), contains 

all the sites of a unique disorder contour in dx , so that the number of disorder contours in dx is 
the same as Ki{S{An+i), X), and we have 

(g^^)-iE-,e8^l7l^Ki(S(A„+i),X) ^ -Q (^75) 

y&dx 

On the other hand, if AT e -^a' +i ' the outermost non-singleton connected component of (S'(A„+i), A") 
has no associated disorder contour in dx, thus 

(^ + ^)-iE^<=a^l7l^«i(S(A„+i),X) -Q ^(^) _ (17g) 

7e9x 



In conclusion, summing over all configurations, we obtain 

^RC.di.ord(A^^^) = (g + ^)^^^^ ^ g-|X|.e(B)-|i3(A„ + ,)\X|.e(0) 



(177) 



XeX^^^"^^ 769 



XeAT^'^'i 769 
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We remark that the definitions (PT|) and may also be extended to general finite volumes of 
the lattice in a compatible fashion. Namely, if for a volume A, we define 

^disord A 1^ c B{A) : dX e Af'°"^} , (179) 
A-^'-d = {x C B{A) : 9 U S(A)^j e A^''^} , (180) 
the compatibility of the definitions can be verified similarly. 

A. 3 A Few Properties of the Biased Random-Cluster Model 

Much information about the standard Potts model can be detected by studying the corresponding 
random-cluster model (see [8j [TOl [H]). Many of the properties of the standard random-cluster 
model can be extended to the biased random-cluster model. These properties, in turn, can be 
used, in a similar fashion, to obtain information about the Potts model with invisible colours. In 
this appendix, we briefly sketch some of these properties that we exploit in the proof of the main 
theorem. The proofs are straightforward modifications of the standard case, which can be found 
in [i [101 [IT]- 

Let G = {S, B) be a finite graph. The configurations of the biased random-cluster model on G 
can be ordered according to the inclusion ordering. A configuration X <Z B is considered to be 
smaller than or equal to a configuration y C B, if and only if every bond present in X is also 
present in Y . An event f C 2^ is increasing if for every two configurations X and Y such that 
X <^ £ and F D A", we have Y <E £. We say that a probability distribution u is positively correlated, 
if nf2) > v{£i)v{£2)- The inclusion ordering on the configuration space 2^ induces an ordering 
on the space of probability distributions on 2^. \i vi and V2 are probability distributions on 2^, 
we write vi < V2 if vi{£) < V2{£) for every increasing event £ C 2^ . In this case, we say that vi is 
stochastically dominated by 1^2- 

For every 0<p<l, q>l and r > 0, the r-biased random-cluster distribution 4>p.q.r on G is 
positively correlated. This follows from the Fortuin-Kasteleyn-Ginibre theorem (Theorem 4.11 of [S]; 
see Corollary 6.7). It follows that, if £ is an increasing (resp., decreasing) event with (j>p^q^r{£) > 0, 
then the conditional distribution (f>p^q^ri' \ £) stochastically dominates (resp., is dominated by) 4>p,q,r- 
Furthermore, if < pi < P2 < I, it follows from Holley's theorem (Theorem 4.8 of [8]) that 
4'pi,q,r ^ 4'p2,q,r (sec CoroUary 6.7 in [8i). 

Let A be a finite volume in the lattice and (jjA the biased random-cluster distribution on A 
(as a graph, without boundary condition). Let us denote by 4>a'^ and the biased random- 

cluster distributions on A with ordered and disordered boundary conditions, respectively. By an 
application of the positive correlation property of 0a we have 

idisord j i j j,ord /-io-i\ 

<Pa '<9k<9k ■ (181) 

Moreover, by a further application of the Fortuin-Kasteleyn-Ginibre theorem, the distributions 
,disord 0ord positively Correlated. This implies that if Ai is a sub- volume of A2, we 



have 



jdisord j idisord j,ord v lord 

CiAi ^ <PA2 and 0^1 ^ 0A2 • (182) 



As in Lemma 6.8 of [H', this implies that the weak limits 



'■'^ ^ lim (/>^''"'^ and ^\iuY(j)°l'^ (183) 

AfL AtL 
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exist, where the hmit A t L can be taken along the net of aU finite volumes in L with the inclusion 
ordering. 

To emphasize the dependence on parameter p, let us write cj)'^'^ and 0a p"^'^ fo^' the biased 
random-cluster distributions with parameter p. Then, by an application of HoUey's theorem, if 
< pi < P2 < 1, we have 

0t;:'^<p;' ^nd rc^^^tp.- (m 

Let i -f^ oo denote the event that there exists an infinite path of bonds passing through site i ( "or- 
der" percolates from site i to infinity). The latter stochastic inequalities imply that the probabilities 
oo) and oo) are increasing in p. This monotonicity assures the existence of 

critical probabilities < < 1 such that for every p < we have (j)™'^{i oo) = 

while for every p > p™^ we have ^ oo) > 0, and similarly for pj?"*"''^. The critical probabil- 

ities are given by 

p"^" ^ sup [p : A (X)) = O} , (185) 

pf^°^'' ^ sup [p : 0f -'^(i A c^) ^ O} . (186) 

It turns out that the two critical probabilities are actually the same, hence we define Pc = p°''^ = 
pdisord^ This follows from the fact that the probability measures (f)'^'^ and 4)^^^°^'^ may differ for at 
most countably many values of p. The latter can be proved in a very similar manner as done in 
Theorem 8.17 of [11] for the standard random-cluster measures. 

By means of the coupling, many properties of the (g, r)-Potts measures can be derived from the 
corresponding r-biased random-cluster measures. For instance, one can show that the thermody- 
namic limits fJ' and f/^^'^ do not depend on the sequence {An}^^ of volumes along which the limits 
are taken. In fact, the limits 

/^lirn^i and i/^^'^ = Ym,^J^ im 

can be taken along the net of all finite volumes in L. In particular, this implies the translation- 
invariance of /i*^ and The proofs are similar to those of the standard case (Proposition 6.9 

of i). 

Uniqueness and multiplicity of the (q, r)-Potts measures are related to the percolation of "order" 
in the r-biased random-cluster model. More specifically, if (f>°'^'^ {order percolate) = 0, then the {q, r)- 
Potts model admits a unique Gibbs measure (as in Theorem 6.10 in ,8 ). On the other hand, if 

/ord f 3 unique infinite sea of order 
\^ with finite islands of disorder 

then the measures /i^, . . . , /i'' are distinct and satisfy 

J. / 3 unique infinite uni-colour sea, 
\^ which has colour k 

(Recall that a "sea" of order in a random-cluster configuration is simply a connected component of 
bonds. A "uni-colour sea" in a Potts configuration refers to a maximal connected subgraph of the 
lattice induced by sites having the same colour.) 



= 1 



(188) 
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The latter claim is a consequence of the existence of a coupling between 0°'''^ and /i'^ (for 
visible fc), which can be constructed as follows: 

i) We first sample a bond configuration X according to 0°''^. 

ii) For every site i that i oo in (S,X), we colour i with colour k. 

iii) For every finite non-singleton connected component of (§,X), we choose a random visible 
colour uniformly among the q possibilities and colour all the sites in the component with this 
colour. 

iv) For every isolated site i in (§, X), we choose a random colour uniformly among the q + r 
possible colours. 

The fact that the marginal of this construction on spin configurations is is parallel to Theo- 
rem 4.91 in [To] and has a similar proof. An analogous coupling exists between and 1^^°'^ . 
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